HU-EP-07/28, DESY 07-108 



The phase structure of a chiraUy invariant lattice Higgs- Yukawa model 

numerical simulations 

P. Gerhold'', K. Jansen'' 
Humboldt-Universitdt zu Berlin, 
Institut fiir Physik, Newtonstr. 15, 
D-12489 Berlin, Germany 
''DESY, 

' Platanenallee 6, D-15738 Zeuthen, Germany 

O ' (Dated: July 25, 2007) 



O 



(N 



in: 



> 

0^ 



The phase diagram ol a chirally invariant lattice Higgs- Yukawa model is explored by means of 
numerical simulations. The results revealing a rich phase structure are compared to analytical large 
Nf calculations which we performed earlier. The analytical and numerical results are in excellent 
agreement at large values of Nf. In the opposite case the large Nf computation still gives a good 
qualitative description of the phase diagram. In particular we find numerical evidence for the 
predicted ferrimagnetic phase at intermediate values of the Yukawa coupling constant and for the 
symmetric phase at strong Yukawa couplings. Emphasis is put on the finite size effects which can 
hide the existence of the latter symmetric phase. 
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I. INTRODUCTION 



The main target of lattice studies of the Higgs- Yukawa sector of the electroweak standard model 
is the non-perturbative determination of lower and upper bounds of the Higgs boson mass |[l|, as 
\ well as its decay properties. There are two main developments which warrant to reconsider these 

OO ■ questions: first, with the advent of the LHC, we are to expect that properties of the standard model 

CO ' Higgs boson, such as the mass and the decay width, will be revealed experimentally. Second, there 

. is, in contrast to the situation of earlier investigations of lattice Higgs- Yukawa models, a consistent 

' formulation of an exact lattice chiral symmetry Q based on the Ginsparg- Wilson relation [4| . 

! Before questions of the Higgs mass bounds and decay properties can be addressed, the phase struc- 

' ture of the model needs to be investigated in order to determine the (bare) couplings in parameter space 

where eventual simulations of phenomenological interest can be performed. There has been a large 
. - activity of investigating lattice Higgs- Yukawa models in the past, see e.g. Refs. [1, S 0, i i, [M [lH 

' for reviews. In particular, the phase structure of lattice Higgs- Yukawa models was investigated in 



. great detail, see e.g. Refs. [12|, |13|, ll4|, [15|, ll6|, |17|, ll8|, |19|, l20| for a still incomplete list. However, in 

these investigations, the lattice formulation of the corresponding Higgs- Yukawa theory broke explicitly 
chiral symmetry. 

This situation changed when it was realized that the Ginsparg- Wilson relation [1] leads to the 
notion of an exact lattice chiral symmetry @ allowing thus to go beyond the earlier models. Based 
on this development, the interest in lattice studies of Higgs- Yukawa models has been renewed [l!], 
H, m, [13, [2^. Here, we follow Liischer's proposition for a chirally invariant and hence consistent 
lattice Higgs- Yukawa model given in Ref. [3]. In particular, we want to address here the question 
whether the phase structure remains as complex as has been found out in the earlier work mentioned 
above. Additional and phenomenologically more interesting questions concerning the behaviour of the 
rcnormalized Higgs and Yukawa couplings will be addressed in future works. 

In 25] we have studied the phase structure of this chirally invariant Higgs- Yukawa model by means 
of a large Nf computation, where Nf denotes the number of fermion generations. We found a com- 
plex phase structure that resembles qualitatively the one of earlier lattice Higgs- Yukawa models, see 
e.g. [14| . which, however, were lacking chiral symmetry. In the present paper we want to confront the 
results obtained in the large Nf approximation with direct numerical simulations for finite values of 
Nf. We remark that in the present work, as in Ref. [1^, we neglect the gauge degrees of freedom and 
consider the pure scalar-fermion sector of the electroweak standard model. 

To be more specific, we consider here a four-dimensional, chirally invariant SU{2)l x SU{2)j^ Higgs- 
Yukawa model discretized on a finite lattice with L lattice sites per dimension such that the total 
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volume becomes V = L^. We set the lattice spacing to one throughout the paper. The model contains 

one four-component, real Higgs field $ and wc consider Nj fermion doublets represented by eight- 
component spinors ■i/i*^*) with i = 1, Nf. Furthermore, there are Nf auxiliary feimionic doublets 
^ serving as a construction tool in the creation of a chirally invariant Yukawa interaction term. 
Once the chiral invariancc is established these unphysical fields arc integrated out leading to a more 
complex model depending only on the Higgs field $ and the Nf physical fermion doublets tp^^\ The 
partition function is written as 

Z= f D^Yl [^"^^'^ ^^^^ ^X^'^ Dx^'^'] exp {-S^ - Sf.'" - Sy) (1) 

i=i 

with the total action being decomposed into the Higgs action S'*, the kinetic fermion action Sp", 
and the Yukawa coupling term Sy- It should be stressed once again that no gauge fields are included 
within this model. 

The kinetic fermion action describes the propagation of the physical fermion fields in the 

usual manner according to 

Nf 

Skin = J2Y1 i^a^T^i^X - 2PX^:^ ln,mXlS (2) 
i=l n,m. 

where the four-dimensional coordinates n, m as well as all field variables and coupling constants are 
given in lattice units throughout this paper. The (doublet) Dirac operator P^™) = P^"") (g) 
is given by the Neuberger overlap operator •£)('"'), which is related to the Wilson operator V^^^ = 
7^^(V^ + V^)-iV^V/by 



= pll + > , i = - p, 1 < p < 2r (3) 




with V^, denoting the forward and backward difference quotients. Note that in absence of gauge 
fields this kinetic part corresponds to the one of free fermions which will be exploited in the numerical 
construction of the overlap operator later. In particular, the eigenvalues of r**^"") can be computed 
analytically. In momentum space with the allowed four- component momenta 

ver-l (-^''^J^' :fori = oo 

^ ^ ^ ~ I {27Tn/L : n e No, n < L}®^ : for i < oo ^ ' 

the eigenvalues of the doublet operator !)('"') are given by 

We remark that the auxiliary fields x^*' do not propagate at all and that their contribution to Sp"' 
is chosen such that the model will obey an exact lattice chiral symmetry. 

The Higgs field couples to the fermions according to the Yukawa coupling term 



(V'^^+X^^) (6) 



n.m i=l 



where yjv denotes the Yukawa coupling constant and will be referred to as Yukawa coupling 

matrix. Here the Higgs field i>„ is rewritten as a quaternionic, 2x2 matrix ^„ = 1 — i^i^Tj, with 
T denoting the vector of Pauli matrices, acting on the SU{2) index of the fermionic doublets. Due to 
the chiral character of this model, left- and right-handed fermions couple differently to the Higgs field, 
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as can be seen from the appearance of the projectors (1 ± 75)/2 in the Yukawa term. Muhiplying 
out the involved Gamma- and Pauh-matrices one can rewrite the couphng matrix in the compactified 
form 



$0l + i$3^5 ^h5+^Kl5 



(7) 



being block diagonal in position space. The model then obeys an exact, but lattice modified, chiral 
symmetry according to 



75 1 



2p 

1 _ J_x)(o^'; 

2p 



(8) 



75 + X^*''75 



recovering the chiral symmetry in the continuum limit [1] . 

Finally, the lattice Higgs action is given by the usual lattice notation 



-K,N 



Nf 



(10) 



n.fi 



with the only particularity that the fermion generation number Nf appears in the quartic coupling 
term which was a convenient convention for the large Nf analysis. However, this version of the lattice 
Higgs action is equivalent to the usual continuum notation 



(11) 



with the bare mass uiq and the bare quartic coupling constant Xq. The connection is established 
through a rescaling of the Higgs field and the involved coupling constants according to 



/K ^ \ 2 1 ^ 2NfXN - 8kn 

(fin = V^KN^n, Ao = — TOn - 



yo 



VN 



where yo denotes the Yukawa coupling constant corresponding to the continuum notation. 



(12) 



II. SIMULATION ALGORITHM 



The first step towards a numerical treatment of the considered Higgs- Yukawa model is to integrate 
out the fermionic degrees of freedom leading to the effective action 

Seff['i>] = S^m- Nf -logdetiM) (13) 

where the fermionic matrix 

M = vnEV^""^ - 2pl?(°'') - 2pyNB (14) 

was given in Ref. [25j . 

Since we focus here on checking the validity of our earlier analytical investigation of the phase 
structure, which was determined in the large iVy-limit, we will only consider even values for Nf, 
allowing to rewrite the effective action according to 

Seff[^ =S^[^~^- logdet {MM"^) , Nf even . (15) 
Thus the positivity of the determinant in Eq. is guaranteed. 
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For the numerical treatment of the remaining determinant in Eq. (|15p we have implemented an 



Hybrid-Monte-Carlo (HMC) algorithm [2€ 
cording to the HMC-Hamiltonian 



|27| |. with Nf/2 complex pseudo-fermionic fields ujj ac- 



Nf/2 



(16) 



where ^ denote the real momenta, conjugate to the Higgs field <&. 

The application of the matrix [A^Aftj-i on ujj can then be performed by means of a Conjugate 
Gradient algorithm due to the hcrmiticity of AiAi^ . However, for the computation of A^x, where x 
is an arbitrary vector, we exploit the fact that there are no gauge fields included within our model. 
The eigenvectors of the used Neuberger overlap operator pf"") are therefore explicitly known to be 
the plane waves 



Xek 



ip), 



Xek 



ip) 



C = ±l, e = ±l, fee {1,2} (17) 



with u'^^{p) denoting the usual four-component spinor structure 




for p ^ and 



for p — Q. 



(18) 



Here G (D^ are two orthonormal vectors and the four component quaternionic vector O is defined 
as G) = (1, it). The corresponding eigenvalues v'^{p) were given in Eq. ([5]). The operators B and 
are thus both block-diagonal, the first in position space and the latter in momentum space. In our 
approach we use a Fast Fourier Transform (EFT) [28] to switch between the position and momentum 
representations, such that all operator applications can be trivially performed due to their actual 
block-diagonal structure. This is particularly advantageous for the overlap operator, since the usual 
construction of this operator would be based on very demanding approximations, e.g. polynomial 
approximations. 

A second advantage of this approach is that the applied Dirac operator can easily be replaced by 
other operators simply by adopting the corresponding eigenvalues. 

Concerning the parallelization of the program there are several options. For example there are 
efficient parallelized EFT-routines available [28|. Here, however, we use a trivial - but very efficient 
- parallelization which is possible due to the large number of fermion generations Nf. We simply 
perform each of the Nf/2 force calculations on a separate computer node. 

For the integration of the obtained forces we find the Leap-Frog integration scheme to be efficient 
on small lattices. This situation changes with increasing lattice size and for L > 16 we get better 
performance with higher order integrators. In that case we use an order 4 Omelyan-integrator [29l.[30j. 
The integration is then performed over a fixed trajectory length set to unity with the typical value 
e = 0.1 for the step size. The step size e is chosen such that the acceptance rate stays between 80% 
and 95%. 

The observables we will be using for exploring the phase structure are the magnetization m and the 
staggered magnetization s, 



m 



i=0 



L4 



■ 3 
El 

.4=0 



4)" 



(19) 



and the corresponding susceptibilities 



Xr. 



V-[{r 



Xs 



(20) 



where (...) denotes the average over the $-field configurations generated in the Monte-Carlo process. 

The auto-correlation of our measurements of these observables in the Monte Carlo time t is then 
accounted for by applying the F-strategy [3l|. In this approach the error a a of an observable A is 
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rewritten as a sum over the correlation function T(t) according to 
r( \ ^ 1 

4 = ciw) = y: m, m = [^^'^ - (^)] • [^^'^'^ ' (^)] (21) 

t=-W ' ' i 

where A^*) denotes the measurement of the observable A in the i-th configuration and N is the total 
number of collected configurations. The variable W is the window in which the function r(t) is to 
be summed up. It should be large enough to obtain reliable estimates of the auto-correlation time r 
which is defined through the exponential decay rate of T{t) 

r(i) ocexp (^-^^ (22) 

and is thus directly connected to the sum C{W) of the auto-correlation function T{t). Typical examples 
for the determination of C{oo) by fitting the function C{W) to a constant are presented in Fig. [TJ 
Since the auto-correlation length r depends strongly on the distance to the phase transition we have 
selected one point in the parameter space close to the phase transition (Fig. [T^) and one point farer 
away from it (Fig. [T]d). Both points correspond, however, to the ferromagnetic phase with a non- 
vanishing Higgs field expectation value, i.e. (to) > 0. For the Higgs field magnetization to as the 

underlying observable we find in these examples the auto-correlation times t™'' = 38.3 ± 1.9 in Fig. [T^ 
and Tm = 7.3 ± 0.4 in Fig.[TjD. We remark here that the value obtained in Fig.[T^ is the largest auto- 
correlation time for the magnetization m encountered in our studies. Although the auto-correlation 
time indeed increases when approaching the phase transition, its value remains acceptable for our 
purposes. Given that our typical statistics is 0(10*) <i>-field configurations, this leads to reliable error 
determinations for the physical quantities of interest. 




FIG. 1: Example for the determination of the auto-correlation time for two different points in the phase 
diagram at Ajv = 0.05, L = 16, Nf = 2, and j/jv = 30. (a): A point very close to the phase transition with 
Kjv = 0.042. (b): A point farer away from the phase transition line with kjv — 0.060. 



III. RESULTS FOR SMALL VALUES OF THE YUKAWA COUPLING CONSTANT 

In this section we will discuss the structure of the phase diagram at small values of the Yukawa 
coupling constant. All numerical results from our simulations are obtained by employing the algorithm 
as detailed in Section|lTl The anticipated structure of the phase diagram can be inferred from our large 
Nf computation in Ref. [1^. In that large Nf approach the Higgs field and the coupling constants 
are scaled according to 

VN = ^= , = ^ , KN = KN , $n = ^/NJ ■ l>„ , (23) 

y/Nf Nf 

where the quantities yN, kn, Aat, and are held constant in the limit Nf — > cxi. Here, we want to 
confront this predicted phase structure with the results of our numerical simulations. 
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At small values of the Yukawa coupling constant, there are two phase transitions when varying kjv: 
the first is a phase transition from a ferromagnetic (FM) phase, with (s) = and (m) > 0, to the 
symmetric (SYM) phase with (m) = (s) =0. The second corresponds to a phase transition from 
the symmetric phase to an anti-ferromagnetic (AFM) phase with (m) — and (s) > 0. These phase 
transitions are expected to be of second order. To locate the phase transition points, we decided to 
fit the data for the susceptibilities Xrm Xs in Eq. (j20p as a function of kjv according to the - partly 
phenomenologically motivated - ansatz 

where A^''^ , ^2*3^1 ^^nd k™;^ are the fitting parameters for the magnetic susceptibility and staggered 
susceptibility, respectively, and 7 denote the critical exponents of the $'*-theory. Here A!^^ {A2 3) is 
actually meant to refer to two parameters, namely for kn < {^n < K^^it) and A™ in 

the other case, such that the resulting curve is not necessarily symmetric. The phase transition point 
is then given at the value of kn = k^^^ {kn = n^rit) where the magnetic (staggered) susceptibility 
develops its maximum. We remark that the ansatz in Eq. (j24p . although not being unique, provides 
a very good description of our numerically obtained data leading to a reliable determination of the 
critical hopping parameters K^j't • 




FIG. 2: An example for the determination of the phase transition points separating the ferromagnetic and 
the anti- ferromagnetic phase from the symmetric phase. We show, as a function of Kiv the behaviour of the 
average magnetization (m) and staggered magnetization (s) in panel (a). The corresponding susceptibilities 
are plotted in panels (b) and (c). The solid lines are fits to the finite size formula of Eq. (|24p. The parameters 
chosen are jjN = 0.632, Ajv = 0.1, L = 6 and Nf = 10. 



In Fig. [2] we present a typical example for the determination of the phase transition points at small 
values of the Yukawa coupling constant. The average magnetizations (m) and (s) as well as the 
corresponding susceptibilities are shown as a function of kjv. We clearly observe the vanishing of the 
magnetization and the staggered magnetization when the symmetric phase is entered (except for some 
small finite volume effects). Associated with these transitions are peaks in the susceptibilities. Note 
that the data for the susceptibilities are fitted very well using the ansatz of Eq. allowing for a 

good determination of the critical points. 



Using the strategy just described we computed the values of k™;^ and n^^.^ for various Yukawa 
coupling constants < 5 holding the quartic coupling = 0.1 constant. In Fig. [3] we summarize 
the numerical results for the phase structure as obtained on 8"*- and 6'*-lattices at Nf = 10 and compare 
them to the analytical Nf = 00, L = 00 phase structure. As expected we observe a symmetric (SYM), 
a ferromagnetic (FM) and an anti-ferromagnetic (AFM) phase, with the symmetric phase bending 
strongly towards smaller values of the critical hopping parameter when the Yukawa coupling constant 
is increased. 
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FIG. 3: The phase diagram at smaU Yukawa couphng constants together with the L = oo prediction of the 
large Nf calculation. The dashed lines denote second order phase transitions while the solid line marks a first 
order transition. The data with open squares were obtained on an 8''-lattice while the ones represented by 
open circles were measured on 6^-lattices. These results were obtained at Ajv = 0.1 and Nf = 10. 

As a general remark we note here that the simulations become extremely demanding when entering 
the anti-ferromagnetic phase, due to an increasingly bad condition number of the fermionic matrix M. 
Within the anti-ferromagnetic phase we thus only present numerical results obtained on 6^-lattices 
throughout this paper. 




-0.3 -0.2 -0.1 0.1 0.2 0.:1 -0.3 -0.2 -0.1 0.1 0.2 0.3 



(a) (b) 

FIG. 4: Evidence for the ferrimagnetic phase with (m) > and (s) > inside the anti-ferromagnetic phase. 
The behaviour of the average magnetization (m) and staggered magnetization (s) is shown in panel (a) as a 
function of hn for yjv = 3.162, Ajv — 0.1, Nf = 10, and L — 6. The corresponding magnetic susceptibility is 
shown in panel (b). From left to right its three observable peaks correspond to the phase transitions AFM-FI, 
FI-AFM, and SYM-FM. From the large Nf, L = oo calculation the ferrimagnetic phase was expected to occur 
approximately at kjv < —0.27. 



Besides these three phases also a fourth, somewhat peculiar phase that can appear at intermediate 
values of the Yukawa coupling constant was predicted by our analytical investigation. This is the so- 
called ferrimagnetic (FI) phase where both, the average magnetization as well as the average staggered 
magnetization, are non-zero, i.e. (m) > and (s) > 0. It was found that such a ferrimagnetic phase 
should exist deeply inside the anti- ferromagnetic phase [l^ . In Fig. 2] we provide evidence for the 
existence of this ferrimagnetic phase. Its location within the phase diagram is in good agreement with 
the analytical prediction. However, the ferrimagnetic phase is not the prime target for our eventual 
interest and hence we do not further investigate this phase here. 

Concerning the order of the encountered phase transitions we find that the SYM-FM as well as the 
SYM- AFM phase transition seem to be of second order in accordance with the continuous behaviour 
of (to) and (s) as seen e.g. in Fig. [4^. This is in contrast to the direct FM-AFM phase transition 
that should occur at intermediate values of the Yukawa coupling constant according to our large Nf 
computation. From the analytical considerations we expect this transition to be of first order. To 
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clarify this we show in Fig. [5] an example for such a phase transition as seen in the numerical simula- 
tions. One can clearly observe an abrupt jump in (m) and (s) in Fig. [5^ indicating a discontinuous 
phase transition. In subfigures (b) and (c) we furthermore present an example for a tunneling event 
between two ground states close to the critical value Kcrit of the hopping parameter serving as another 
strong indication for the first order nature of the phase transition at intermediate values of the Yukawa 
coupling constant. However, we do not study the order of the phase transition in great detail here, 
since this is not in our main interest. 




-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 (j 2OOO 4000 0000 8000 lOOOO 2000 4000 0000 8000 10000 

f^jV Monte Carlo Time Monte Carlo Tinre 

(a) (b) (c) 

FIG. 5: The direct FM-AFM phase transition at intermediate values of the Yukawa coupling constant. We 
show, as a function of kjv, the behaviour of the average magnetization (m) and staggered magnetization (s) 
in panel (a). The parameters chosen are yjv = 6.325, Ajv = 0.1, L = 4, and Nf — 10. Panels (b) and (c) show 
a tunneling event between two ground states which we take as strong indication for the first order character of 
the phase transition. The plots show m in panel (b) and s in panel (c), respectively, versus the Monte-Carlo 
time at the hopping parameter kjv ~ —0.196 being very close to its critical value. 



Qualitatively, all presented findings are in excellent accordance with our large Nf calculations in 
Ref. 25]. On a quantitative level, however, the encountered deviations in Fig. [3] need to be further 
addressed. These deviations can be ascribed to finite volume effects as well as finite Nf corrections. 
Here we start with a discussion of the finite volume effects. 

The location of the phase transition points can be strongly altered by finite size effects. This is 
illustrated in Fig. [5] showing some phase transition points from the FM to the SYM phase as obtained 
from our numerical simulations on a 4^-lattice (open squares), and on an B^'-lattice (open circles). 
One clearly observes that the phase transition line is shifted towards smaller values of the hopping 
parameter when the lattice size is increased. This effect can also be anticipated from the analytical 
computation of the phase transition line, when one imposes finite lattices also for the minimization 
of the effective potential in the large Nf approximation. Since we want to demonstrate the finite 
volume dependence here isolated from the iV/-dependence, we present the numerical results for the 
(very large) value of fermion generations Nf = 50 and compare them to the analytical Nf = oo phase 
transition lines obtained for L = 4 (dotted line), L = 8 (dashed line), and L = oo (solid line). The 
analytical lines perfectly describe the numerical results and one clearly observes the convergence of 
the numerical results to the analytically predicted L = oo line as the lattice size increases. 

However, one remark is in order here for the orientation of the reader, which concerns the large 
TVj^-computation of the phase transition points in a finite volume: The fermionic determinant det(M) 
with M given in Eq. p4|) becomes, on a finite lattice, identical to zero for completely vanishing Higgs 
field. On infinite lattices the zero modes of V^"""^ form a set of only zero measure and the integral 
entering the effective action can be shown to converge, such that there actually is a symmetric phase 
on infinite lattices. For finite L we therefore cannot determine the phase transition by simply searching 
that value of the hopping parameter, where the average magnetization vanishes. Instead we search 
for that K^v, where the minimum of the effective action Seff becomes flattest, i.e. where the second 
derivative of Sef / with respect to the magnetization becomes minimal at the location of the minimum. 
Since the Higgs field oscillates the stronger around the minimum of the effective action the smaller its 
second derivative is, this approach corresponds to finding the phase transition point by searching for 
the maximum of the susceptibility. 
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FIG. 6: A demonstration of finite size effects. We show for several selected values of the Yukawa coupling 
constant the phase transition points between the ferromagnetic and the symmetric phase as obtained on a 
4*-Iattice (open squares), and on an 8*-lattice (open circles). These results are compared to the analytical 
1/ = 4 (dotted), L = 8 (dashed), and L = oo (solid) phase transition lines determined in the large A^'^-limit. 
The chosen parameters are Ajv = 0.1 and Nf = 50. 



The A'^/-dependence of the numerically obtained critical hopping parameters and K^^jj is shown 
in Fig. [7| for several selected values of the Yukawa coupling constant. One clearly sees that for 
increasing Nf the numerical results converge very well to the analytical finite volume predictions, as 
expected. It is interesting to note that the leading term in the finite Nf corrections, i.e. the 1/Nf 
contribution, seems to be the only relevant correction here, even at the small value Nf = 2, as can be 
seen in Fig. [7] by fitting the deviations to the function fm,siNf) = Am,s/Nf with Am.s being the only 
free parameter. Furthermore, one observes that the critical hopping parameter k™;^. is shifted towards 
larger values with decreasing Nf while is shifted towards smaller values. 
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FIG. 7: The Af/-dependence of the critical hopping parameters K^it) '^crit for the selected Yukawa coupling 
parameters j/jv = 0.0 (a), yjv = 1.0 (b), and j/jv = 2.0 (c). The data with square symbols were measured on 
an 8*-lattice while those represented by circles were obtained on 6*-lattices. The analytical, finite volume, 
large Nj predictions for the SYM-FM (SYM-AFM) phase transitions are represented by the dashed (dotted) 
lines. The dash-dotted lines are fits of the numerical data to the function fm,s{Nf) — Am.s/Nf + Bm,s where 
Bm,3 is set to the actual analytical prediction and Am.s is the only free fitting parameter. The results were 
computed for A;v = 0.1. 



From our findings in this section we finally conclude that the structure of the phase diagram of 
the considered Higgs- Yukawa model at small values of the Yukawa coupling constant can be very well 
predicted on a qualitative level by the results of our large Nf analysis. It also gives a very good 
understanding of the encountered finite volume effects. 
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IV. RESULTS FOR LARGE YUKAWA COUPLING CONSTANT 



In this section we want to address the region of large Yukawa coupHng constants, i.e. j/a? ^ 1. From 
our large Nf calculations we expect here a ferromagnetic, an anti-ferromagnetic and a symmetric 
phase. The large Nf calculation also revealed that significant finite size effects can be present in the 
symmetric phase which may render its detection difficult. This large Nf approach was carried out by 
scaling the Higgs field and the coupling constants according to 

VN =yN, Aat = ^ , KAT = , $„ = y/Nj ■ l>„ , (25) 

Nf Nf 

where the quantities y^, A^r, kjv, and were held constant in the limit Nf oo. 




FIG. 8: The behaviour of the average magnetization (m) and staggered magnetization (s) as a function of kn 
on a 4''- (a), S**- (b) and 16'*-lattice (c). In the plots we have chosen yjv = 30, Ajv — 0.1 and Nf = 2. 



In Fig. [51 we show the numerically obtained values for the average magnetizations (m) and (s) on 
various sized lattices as a function of kn for a large value of the Yukawa coupling constant j/at — 30. 
Fig. [5] demonstrates that indeed the symmetric phase emerges only on sufficiently large lattices, while 
on small lattices the magnetization does not vanish as a function of decreasing kn even deeply within 
the anti-ferromagnetic phase. Instead (m) reaches a plateau with a clearly non-vanishing value in the 
limit Kjv — > — oo. This becomes especially well observable for the smallest considered lattice, the 4*- 
lattice presented in Fig. [5^. Thus, one may erroneously conclude that there is no symmetric phase at 
large values of the Yukawa coupling constant, if one considers too small lattices. However, the plateau 
value of (m) is fully consistent with our analytical results predicting a finite volume effect causing a 
non- vanishing magnetization (m) > also for arbitrarily negative values oi kn- To demonstrate this 
latter statement we restate here one result of Ref. ,25] for the effective action of a field configuration 
in terms of its magnetizations m and s in the large ?/Ar-limit, reading 

5e//[$] =5*-iV/-^81og 

n 

with the abbreviations 

rh = — 5 ^ and s — — -. (27) 

— s — m'' 

Considering only the ground state of this effective action one cannot correctly predict the phase 
transition of the model, as discussed in Ref. [l^. However, it is sufficient to correctly predict the 
behaviour of (to) and (s) in the limit of large negative (and positive) values of the hopping parameter 
Kn, as demonstrated in Fig. [9l where we plot again the average magnetizations for the 4'^-lattice 
together with the finite volume analytical expectations, obtained by minimizing the effective action of 
Eq. p6p. The convergence of the numerical results to the analytical finite volume prediction is very 
well observed in Fig. 



+ s.(-l) 



7V/-81og|m|-iV/-561og 



Ito^ 



(26) 
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FIG. 9: Comparison with finite size expectations. We show again the average magnetizations (m) and (s) of 
Fig. UK a-s obtained on the 4^-iattice at jji^ = 30, Ajv = 0.1, and Nf — 2. We compare these magnetizations to 
the analytical predictions for (m) (dashed line) and (s) (dotted line) computed by minimizing the effective, 
large j/jv action of Eq. (|26|l . Panel (b) is just a magnification of plot (a). 



We remark that the non- vanishing plateau as well as the asymmetry in (m) and (s) are both caused 
by the term log \rh\ appearing in Eq. (I26p . This term as well as the very last one in this equation do, 
however, not scale proportional to the volume L* in contrast to all other contributions to the effective 
action. Its influence therefore eventually disappears as the lattice size increases. This is exactly what 
is observed here. 




KJV f-'A- 

(a) (b) (c) 

FIG. 10: The behaviour of the magnetic susceptibility Xm as a function of the hopping parameter kjv on a 
4*- (a), S**- (b) and 16^-lattice (c). In the plots we have chosen yjv = 30, Ajv — 0.1 and Nf — 2. The fit in 
panel (b) is only applied to those points with kjv > 0.025 or kjv < —0.05 in order to reduce the influence of 
the unphysical peak at kjv = 0.0. Note the changing scale in the three plots. 



In Fig. [To] we show the susceptibilities Xm corresponding to the magnetizations in Fig. [S] For the 
smallest lattice, i.e. the 4*-lattice, one observes only one peak in the magnetic susceptibility, centered 
at kn = 0. From this result one could conclude that the phase transition point is located at kjv = 0, 
excluding a symmetric phase, since the staggered susceptibility reaches its maximum at the same value 
of kat. However, with increasing lattice sizes a second peak develops in the susceptibilities. This is 
very well observed in Fig. corresponding to the larger 8'*-lattice. It shows that indeed two distinct 
peaks emerge on this intermediate lattice. It is actually this second peak, centered around kn = 0.04 
in this case, that correctly describes the physical phase transition between the ferromagnetic and the 
symmetric phase, while the first one is only caused by the finite volume terms discussed in |25| . which 
do not scale with the lattice volume. Its height is therefore at most constant in contrast to the physical 
peak, which grows with increasing lattice volume. On the largest presented lattice, the 16'*-lattice, 
the physical peak a,t kn = 0.04 completely dominates the scene and the former small volume peak at 
kn = has disappeared, presumably hidden beneath the large error bars at kjv — 0. 
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FIG. 11: The phase diagram for large values of the Yukawa coupling constant. We show the numerical results 
for the phase transition points to the ferromagnetic phase as obtained on an 8*-lattice at Ajv = 0.1, and 
Nf = 2. We could not reliably determine the phase transition points to the anti-ferromagnetic phase, as 
explained in the main text. We compare these results to the prediction of our large Nf calculation for the 
SYM-FM transition (dashed line). The SYM- AFM phase transition is marked by the dotted line. 

We have then determined the SYM-FM phase transition points by fitting the physical peaks in the 
magnetic susceptibility Xm on the intermediate 8'*-lattices to the finite volume expectation in Eq. (I24p 
by taking only the points belonging to the physical peak into account as demonstrated in Fig. llOb . 
We remark here that we do not provide any data for the SYM-AFM phase transition because the 
phase transition is not reliably detectable on the 6'*-lattice, due to the finite volume effects discussed 
above, and the 8*-siniulations are not practicable in the anti-ferromagnetic phase with our algorithm 
as explained in Section Hill 

In Fig. [11] we finally summarize the obtained phase transition points together with the analytical 
Nf — 00, L — 00 expectation of the phase structure at large values of the Yukawa coupling constant. 
Qualitatively, the picture we obtain from the numerical simulations is in full accordance with the results 
from the large Nf approximation: there are second order phase transitions separating a ferromagnetic 
phase from a symmetric phase. In this symmetric phase strong finite size effects are encountered such 
that only for large lattice sizes this symmetric phase can be identified. Quantitatively, the numerical 
results deviate from the analytical expectation due to the finite settings Nf = 2 and L = 8, but are 
still in good agreement. 

In this section we found that a symmetric phase at large values of the Yukawa coupling constant 
does indeed exist although its existence is obscured on too small lattices by strong finite size effects, 
and that its location within the phase diagram is in good agreement with the analytical large Nf 
predictions. We remark here that the existence of a symmetric phase at strong Yukawa couplings has 
also been observed and discussed in Ref . [l^l • From our findings we thus conclude that the analytical 
large Nf calculations describe the phase structure of the considered Higgs- Yukawa model at large 
values of the Yukawa coupling constant very well. 

V. SUMMARY AND OUTLOOK 

In this paper we have studied by numerical simulations the phase structure of a chirally invariant 
lattice Higgs- Yukawa model, originally proposed by Liischcr, in order to check the validity of our 
earlier analytical investigation of its phase structure. These earlier calculations have been performed 
in the large A^/-limit for small and for large values of the Yukawa coupling constant. 

In Section IIIII we compared the numerical to the analytical results at small values of the Yukawa 
coupling constant. We started with a discussion of the qualitative structure of the phase diagram. For 
that purpose we presented our numerical results for the phase transition lines obtained at Nf = 10 on 
some 8^- and 6^-lattices and compared them to the analytically computed N f ~ 00, L = 00 predictions. 
Qualitatively, the numerical and analytical results are in very good agreement: As expected we clearly 
observe a symmetric (SYM), a ferromagnetic (FM) and an anti-ferromagnetic (AFM) phase. With 
increasing Yukawa coupling constant the symmetric phase strongly bends downwards to smaller values 
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of the hopping parameter. In particular, we find that the obtained phase structure resembles that of 
earlier Higgs- Yukawa models on a qualitative level. One peculiarity in these types of models is the 
emergence of a ferrimagnetic (FI) phase with (m) > and (s) > deeply inside the anti-ferromagnetic 
phase. This phase was also predicted by the analytical investigation of the model and it is located 
at the predicted position within the phase diagram. Furthermore, the data also support very well 
the analytical expectations concerning the order of the encountered phase transitions. The SYM-FM 
and the SYM-AFM phase transitions were supposed to be of second order while the direct FM-AFM 
transition was predicted to be of first order. Although we did not study that in great detail, since 
this aspect was not in out main interest, the obtained lattice results are in very good agreement with 
these analytical findings. 

We then turned towards the quantitative discussion of the encountered deviations between the 
numerical finite volume, finite Nf results and the presented analytical Nf = oo, L — oo calculations. 
Firstly, we showed that finite volume effects alter the location of the phase transition lines strongly. In 
order to isolate the finite volume effects from the 1/iV/ corrections we presented numerical results for 
the phase transition points for the choice of the very large number of fermion generations Nf = 50. In 
that setting we could show that the finite volume effects are in excellent agreement with the analytical 
finite volume predictions. We then demonstrated the strength of the 1/Nf corrections by presenting 
the numerically obtained phase transition points at smaller values of Nf for some selected Yukawa 
coupling parameters. We found that the 1/-/V/ corrections drive the critical hopping parameters 
towards larger values for the case of the SYM-FM phase transition and towards smaller values for the 
SYM-AFM transition. Besides that these corrections do not change the qualitative phase structure of 
the model. 

We then discussed the phase structure at large values of the Yukawa coupling constant in Section HVl 
In particular we showed that there is actually a symmetric phase in this regime of the Yukawa cou- 
pling constant and that it is located at the expected position within the phase diagram. We also 
demonstrated that this symmetric phase becomes unobservable on too small lattices due to strong 
finite volume effects, as derived in our earlier studies, preventing the Higgs field expectation value 
from vanishing. We furthermore showed that the behaviour of the magnetization at large negative 
(and positive) values of the hopping parameter km can be very well described by taking these finite 
volume contributions into account. The emergence of the symmetric phase with increasing lattice 
size could also clearly be observed in the presented plots of the magnetic susceptibility Xm (Fig. [TU|) . 
Finally, we presented our numerical results for the critical hopping parameters k;™j^ of the SYM-FM 
phase transition at large Yukawa coupling constants. We compared them to the analytical large Nf 
predictions and found them to be in good agreement even though the numerical simulations were 
performed at Nf = 2. 

We end with a short outlook about our next steps concerning the further investigation of the 
presented Higgs- Yukawa model: We have started the implementation of a PHMC algorithm [s^ with 
which the simulation will become possible at arbitrary values of Nf, in particular at the physically 
interesting setting Nf = 1. Having the qualitative phase diagram of the model at hand we will then 
search for the physical region of the parameter space, reproducing the top quark mass, eventually 
allowing to find upper and lower bounds for the Higgs boson mass. 
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